Dynamic Alignment in Driven Magnetohydrodynamic Turbulence 
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Motivated by recent analytic predictions, we report numerical evidence showing that in driven 
incompressible magnetohydrodynamic turbulence the magnetic- and velocity-field fluctuations lo- 
cally tend to align the directions of their polarizations. This dynamic alignment is stronger at 
smaller scales with the angular mismatch between the polarizations decreasing with the scale A 
approximately as 9x oc A^''*. This can naturally lead to a weakening of the nonlinear interactions 
and provide an explanation for the energy spectrum E(k) oc k'^^^ that is observed in numerical 
experiments of strongly magnetized turbulence. 

PACS numbers: 52.30.Cv, 52.35.Ra, 95.30.Qd 



INTRODUCTION. 

Incompressible hydrodynamic turbulence is described 
by the Kolmogorov theory of energy cascade, in which 
turbulent fluctuations (or eddies) transfer energy from 
the largest scales, where it is injected, to smaller ed- 
dies until the dissipative scale is reached and energy is 
extracted from the system (e.g., [H [I])- Such an en- 
ergy cascade is local, i.e. its rate at each scale depends 
solely upon the characteristics of eddies at that scale. 
It follows that an eddy of size A, with typical velocity 
fluctuations of strength Sv\, will lose energy to (or get 
fragmented into) smaller-scale eddies over a time of du- 
ration T\ oc X/dv\ (the only quantity having the dimen- 
sions of time that can be constructed in this system). In 
a steady state, the energy flux from the injection scale 
to the dissipative scale is the same at all scales, leading 
to 5v\/tx = const and 5vx cx A^/"^. The corresponding 
Fourier spectrum of velocity fluctuations then has the 
Kolmogorov scaling, Ek cx k^^^^ . 

Iroshnikov 0] and Kraichnan [J| realized that magne- 
tohydrodynamic (MHD) turbulence is qualitatively dif- 
ferent from non-magnetized turbulence. The governing 
equations of incompressible MHD are (see, e.g., Q) 

(9tv -I- (v • V)v = -Vp + (V X B) X B + i/Av + f , 

5tB = V X (v X B) + 77AB, (1) 
V-v = 0, V-B = 0, 

where v(x, i) is the velocity field, B(x, t) the magnetic 
field, p the pressure and f (x, t) is the external force. Un- 
like the large-scale velocity which can be removed from 
the hydrodynamic system by means of a Galilean trans- 
formation (thus allowing the hydrodynamic system to be 
treated locally), here the magnetic field of large-scale ed- 
dies cannot be eliminated by transforming into a moving 
reference frame. Thus the small-scale eddies always ex- 
perience the action of the large-scale magnetic field. The 
MHD turbulent cascade is therefore mediated by such a 
guiding field. 



The Iroshnikov-Kraichnan theory of MHD turbulence 
is formulated upon the assumption that the turbulence is 
three-dimensionally isotropic. In this theory, the turbu- 
lence becomes progressively weaker as the cascade pro- 
ceeds to smaller scales and the turbulent spectrum has 
the form Ek oc A detailed discussion of the 

Iroshnikov-Kraichnan theory can be found, for example, 
in [2]. ^ 

The isotropy of turbulence in the presence of a strong, 
large-scale magnetic field was questioned in early ana- 
lytic and numerical considerations (see, e.g., @, [1,13, [1] ) , 
and it was demonstrated that the energy cascade is di- 
rected mostly perpendicularly to the guiding field (e.g., 
P]). In 1995, Goldreich and Sridhar [l^] developed a 
theory by proposing that the turbulent eddies become 
progressively more stretched along the guiding field as 
the cascade proceeds toward small scales. As a result, 
the time of nonlinear interaction can be estimated as in 
the Kolmogorov theory, however, the energy cascade is 
anisotropic with respect to the guiding field. The field- 
perpendicular energy spectrum then has then the form 
E{ki) oc k'^'"^. 

In recent years, increasingly high-resolution numeri- 
cal simulations have indeed confirmed the anisotropy of 
MHD turbulence. The turbulent fluctuations are elon- 
gated along the guiding field (see, e.g., 11, l^, 13. l3|). 
and the anisotropy is scale-dependent. However, the nu- 
merical simulations also find the field-perpendicular en- 
ergy spectrum E{kA_) oc k'^'"^ ([IlQIli; see also [11]). 
Thus the findings are neither described by the isotropic 
Iroshnikov-Kraichnan theory, nor do they agree with the 
Goldrcich-Sridhar scaling for the energy spectrum. 

A possible resolution to this controversy has been re- 



cently proposed in [17|, |18|. To discuss it, let us first 
note that at each wave number k, the Fourier compo- 
nents of fluctuating fields can be expanded in shear- 
Alfven and pseudo-Alfven modes (we will provide more 
detail in the next section). Since the turbulent cascade 
proceeds mainly in the field perpendicular direction, it 
is dominated by shear-Alfven modes, while the pseudo- 
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Alfven modes are passively advected by the turbulent 
cascade Cs ee lioll fo r a detailed discussion). It was sug- 
gested in [iTI. Il8l| that the polarizations of the shear- 
Alfven magnetic-field and velocity-field fluctuations be- 
come spontaneously aligned in a turbulent cascade. This 
alignment becomes progressively stronger for smaller 
scales. It was proposed that at a given scale A the typical 
fluctuations dvx and ±Sh\ are aligned within the angle 



i\ oc 



(2) 



As in the Goldreich-Sridhar theory, the eddies are 
stretched along the guiding field. However, as a result 
of polarization alignment, they also become anisotropic 
in the field-perpendicular plane, and this anisotropy in- 
creases as the scale decreases. This leads to a scale- 
dependent depletion of the nonlinear interaction in ([T]) 
and to the spectrum of field-perpendicular fluctuations 
E{k±) cx kj_^^^, in good agreement with the numerical 
results described above. 

As proposed in 17 , 11] , the reason for such a dynamic 



alignment may be the existence of two conserved quan- 
tities in magnetohydrodynamics, whose cascades are di- 
rected toward small scales in a turbulent state. The MHD 
system ([T]) in the ideal limit ly — rj = conserves the in- 
tegral of energy 
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(6^ + v^)£x, 



and the integral of cross-helicity, 



^ y"(v • h)d^x, 



(3) 



(4) 



provided that the fluctuations v(x) and b(x) have pe- 
riodic boundary conditions or vanish at inflnity. It was 



proposed in [17|, [18| that the sole requirement of con- 



stant energy flux does not allow one to flnd the spectrum 
of the fluctuations uniquely. The structure of turbulent 
fluctuations in the inertial region should accommodate 
constant fluxes of both conserved quantities, which leads 
to scale-dependent anisotropy of turbulent eddies in the 
field perpendicular plane, and to the alignment 

Interestingly, the effect of dynamic alignment has been 
extensively investigated in the case of decaying MHD 
turbulence (e.g. [ll, [H, |2l|, [23, HI]), where it essen- 
tially means that decaying magnetic and velocity fields 
approach asymptotically in time the so-called Alfvenic 
state v(a;) ~ ±b(a;). However, decaying turbulence is 
qualitatively different from its forced counterpart, and 
the effects that we discuss in the present work have not 
been addressed in the earlier investigations. 

We also mention that previous explanations of the 
numerically observed shallower-than-Kolmogorov spec- 
trum have essentially invoked intermittency effects, e.g., 
[3 Hi ESI- Although intermittency may significantly 
affect the scaling of higher-order correlation functions. 



it usually provides only small corrections to the energy 
spectrum. Our explanation of the -3/2 spectrum does 
not require intermittency. 

In the present paper we investigate the phenomenon of 
dynamic alignment via direct numerical simulations. We 
analyze driven incompressible MHD turbulence with a 
strong guiding magnetic field. We measure the degree to 
which the velocity and magnetic fleld fluctuations align as 
a function of scale and we also investigate the dependence 
on the strength of the imposed fleld. Numerical veriflca- 
tion of the scale-dependent dynamic alignment Q is the 
main goal of our work. 



DYNAMIC ALIGNMENT IN MHD 
TURBULENCE. 

We solve the MHD equations ([T]) using standard pseu- 
dospectral methods. An external magnetic fleld is ap- 
plied in z direction with strength Bq measured in units 
of velocity. The periodic domain has a resolution of 256"^ 
mesh points and is elongated in the z direction, with as- 
pect ratio l:l:i?o- The external force, f (x, t), is chosen so 
as to drive the turbulence at large scales and it satisfles 
the following requirements: it has no component along z, 
it is solenoidal in the x — y plane, all the Fourier coeffl- 
cients outside the range 1 < fc < 2 are zero, the Fourier 
coefficients inside that range are Gaussian random num- 
bers with amplitude chosen so that the resulting rms ve- 
locity fluctuations are of order unity, and the individual 
random values are refreshed independently on average ev- 
ery turnover time of the large scale eddies. The Reynolds 
number is defined as Re — UrmsL/v, where L 1) is 
the fleld-perpendicular box size, v is fluid viscosity, and 
Urms 1) is the rms value of velocity fluctuations. We 
restrict ourselves to the case in which magnetic resistiv- 
ity and fluid viscosity are the same, v = rj. The system 
is evolved until a stationary state is reached (conflrmed 
by observing the time evolution of the total energy of 
fluctuations) and the data are then sampled in intervals 
of order of the eddy turn-over time. All results presented 
correspond to averages over these samples (approx. 10 
samples). 

First we measure the two-dimensional energy spec- 
trum, deflned as E{ki_) = {\w{k±_)\'^)k±_ + {\h{k±_)f)ki^, 
where v(fcj^) and b(fcx) are two-dimensional Fourier 
transformations of the velocity and magnetic flelds in a 

1/2 

plane perpendicular to Bq and k± — (fc^ + fc^) . The 
average is taken over all such planes in the data cube, 
and then over all data cubes. The resulting spectrum of 
fluctuations is presented in Fig. [TJ It is impossible to in- 
fer the exponent of the power-law distribution with good 
accuracy here and in particular it is hard to distinguish 
between the spectral indices —5/3 and —3/2. As was 
demonstrated in a much higher resolution al- 

lows greater Reynolds numbers to be explored and yields 
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FIG. 1: Two-dimensional spectrum of MHD turbulence in the 
plane perpendicular to the guiding field Bq. The numerical 
resolution is 256"^, the Reynolds number is Re = 800, and 
the strength of the uniform guiding field is Bo ~ 5. The 
turbulence is driven by a solenoidal random force at wave 
numbers = 1, 2, as is explained in the text. 



a longer inertial range. Since recent high resolution cal- 
culations do yield E{k±) - kj^^^ (see 0,111 EH) we 
choose not to pursue this issue further here. Instead we 
concentrate our study on the effect of dynamic alignment 
that, as we shall demonstrate presently, can be observed 
well even with limited resolution. 

To investigate the dynamic alignment of the shear- 
Alfven fluctuations in the field-perpendicular plane, we 
use specially constructed structure functions. Let us de- 
note 6vr = v(x + r) — v(x) and 6hr = b(x + r) — b(x), 
where r is a point-separation vector in a plane perpendic- 
ular to the large-scale field, Bq (see [1^). In incompress- 
ible MHD there are two types of linear modes: shear- 
Alfven waves and pseudo-Alfven waves. Consider a wave 
whose wave vector k is almost perpendicular to the guid- 
ing field. Then the shear- Alfven wave will have polariza- 
tion of V and b fluctuations in the direction perpendicular 
to both Bq and k. The polarization of the pseudo-Alfven 
wave will lie in the plane of Bq and k, perpendicularly 
to k. Since the wave vector is almost perpendicular to 
Bq, the polarization of the pseudo-Alfven wave will be 
closely aligned with Bo. 

In the nonlinear case, the Fourier component of fluc- 
tuations at each wave vector k can be expanded into 
two components: those with polarizations in the Bq x k 
direction, and those in the (Bq x k) x k direction. Al- 
though, in contrast with the linear case, such modes are 
generally not solutions of the MHD equations, we will 
refer to them as the shear- Alfven mode and the pseudo- 
Alfven mode, respectively. As we mentioned earlier (see 
also the shear- Alfven modes play the dominant 

role in the anisotropic cascade, while the pseudo-Alfven 
modes are passively advected by them. The scaling of 
the pseudo-Alfven spectrum then follows the scaling of 



the shear- Alfven spectrum. However, since the pseudo- 
Alfven modes do not dominate the dynamics, they should 
be excluded when the angular alignment is calculated. In 
other words, in order to restrict ourselves to shear- Alfven 
fluctuations we need to exclude the components of (5v,. 
and 5hr parallel to the large scale magnetic field. 

It is important to note however, that since the tur- 
bulence is strong an eddy of size r lives only one eddy 
turn-over, during which time it can be transported along 
the large scale field only by a distance comparable to its 
own size. Therefore, during its life time, such an eddy 
feels only the local direction of the guiding field which can 
be tilted with respect to the direction of the global field 
Bq by the action of larger and longer lived eddies; this 
was recognized by Cho & Vishniac ll| . Therefore, when 
calculating the angular alignment we need to remove 
the components of Sv^ and Shr that are in the direc- 
tion of the local magnetic field B(x) (when Bg is strong, 
B(x) « Bq). Thus we calculate Svr — (5v,. — {Sv,. ■ n)n 
and 5hr — Shr — {6hr ■ n)n, where n = B(x)/|B(x)|. 

We are now ready to introduce the following two struc- 
ture functions. The first is defined by 



,{r) = {\Svr X 5hr\), 



(5) 



where x denotes a vector cross-product, and the average 
is performed over all different positions of the point x. 
The second structure function is defined by 



S^ir) = (|5v,.||<5b,|). 



(6) 



with the same averaging procedure. By definition of the 
cross product, the two structure functions are related by 
the angle between the polarizations of and (5b,. . This 
is precisely the phenomenon we would like to investigate. 
When the alignment angle Or, say, is small, we have 



sm ( 



,.) = Scrossir)/S2{r). 



(7) 



Presented in Fig. [5] is the result of such a numerical cal- 
culation. There the angle is compensated by the phe- 
nomenological scaling r^^^ of 13, T^- The data shows 
a clear good agreement with the theory from approxi- 
mately 1/50 to 1/10 of the field-perpendicular box size. 
The effects of varying the guiding field strength are shown 
in the insert of Fig. [2l We note that provided that the 
field is sufficiently strong the alignment effect is robust 
and does not change significantly with field strength. For 
weaker fields on the other hand, the alignment is gradu- 
ally lost. 

We should note that recently, Beresnyak and Lazar- 
ian [2J] atteiripted to measure the geometric alignment 
predicted in [U [l3| by a procedure analogous to av- 
eraging the angle between the normalized fluctuations 
Svr = S-Vr/\S-Vr\ and Shr = Shr/\Shr\. They did not 
observe the alignment ^ and concluded that the weak- 
ening of interaction was, rather, an intermittency effect. 



4 



CD 



1.3 

LO 
0.9 

0.8 
0.7 









0.30 


















0,25 


+ + 








0.20 


+ 






f/; 


0.15 


f 








0.10 










2 4 6 8 10 12 

1 = 



0.01 



0.10 



FIG. 2: The alignment angle, defined in Eq. ((Tj, compensated 
by the theoretical scaling r^l^ {Bo = 5, i?e = 800). The 
insert demonstrates how the slope d9/dr in the inertial range 
changes with the strength of the guiding field. 



Their averaging procedure, however, does not respect the 
fact that, by definition, the 5-vx and 5hx to which 
refers are typical fluctuations at the scale A, that is, fluc- 
tuations whose strengths are close to their rms values (cf. 
p., 2]). The alignment angle should be measured between 
such dynamically relevant fluctuations, as, for instance, 
in our approach ([7|). 

The numerical verification of the scale dependent dy- 
namic alignment of the magnetic and velocity polariza- 
tions in driven MHD turbulence, presented in Fig. [2l is 
the main result of our work. 



DISCUSSION 

Two consequences of the observed polarization align- 
ment have particular importance. These concern the en- 
ergy spectrum and the viscous scale (or inner scale) of 
MHD turbulence. 

Energy spectrum. — The scale dependent alignment can 
naturally imply that the energy transfer time is t\ 



(see 



171 , Il8l|). Since we obtained 9\ oc 
A^/*, the requirement of constant energy flux, 5v\/tx = 
const, then leads to 5vx oc A^/*, where A is the field- 
perpendicular scale of fluctuations. This translates to 
the field-perpendicular energy spectrum of MHD turbu- 
lence, E{k±) oc kj^^"^ , announced in the introduction. 

Viscous scale. — At the viscous scale Aj/, the time of 
nonlinear interaction tx is of the order of the diffusive 
time Tj/ ~ }? jv. A simple calculation then leads to \y ^ 
L/Re^^"^. This result is qualitatively different from that 
for nonmagnetized turbulence, Ay(i? = 0) ~ L/Ke^^^ 
(see e.g., [l|). This implies that for the same Reynolds 
numbers, the viscous scale in MHD turbulence is larger 
than that in nonmagnetized turbulence. This difference 



is especially relevant for astrophysical plasmas, where 
Reynolds numbers are quite large. Re ~ 10^ — 10^". 

Apart from its fundamental value, the existence of dy- 
namic alignment in driven MHD turbulence has conse- 
quences for our understanding of such astrophysical phe- 
nomena as solar wind structure, interstellar scintillation, 
cosmic-ray transport in galaxies, and heat conduction in 
galaxy clusters. A discussion of these matters will be 
presented elsewhere. 
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